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In spin values the damping parameters of the free layer are determined non-locally by the entire 
magnetic configuration. In a dual spin valve structure that comprises a free layer embedded be- 
tween two pinned layers, the spin pumping mechanism, in combination with the angular momentum 
conservation, renders the tensor-like damping parameters tunable by varying the interfacial and dif- 
fusive properties. Simulations based on the Landau-Lifshitz-Gilbert phenomenology for a macrospin 
model are performed with the tensor-like damping and the relaxation time of the free layer mag- 
netization is found to be largely dependent on while tunable through the magnetic configuration of 
the source-drain magnetization. 

PACS numbers: 75.70.Ak, 72.25.Ba, 75.60.Jk, 72.25.Rb 



A thorough knowledge of magnetization relaxation 
holds the key to understand magnetization dynamics in 
response to applied fields^ and spin-transfer torquesj^ 
In the framework of Landau-Lifshitz-Gilbert (LLG) phe- 
nomenology, relaxation is well captured by the Gilbert 
damping parameter that is usually cited as a scalar quan- 
tity. As pointed out by Brown half a century agoj^ the 
Gilbert damping for a single domain magnetic particle is 
in general a tensor. 

When a ferromagnetic thin film is deposited on a nor- 
mal metal substrate, an enhanced damping has been ob- 
served ferromagnetic resonance experiments.^ This ob- 
servation is successfully explained by spin pumpingi^i 
The slow precession of the magnetization pumps spin 
current into the adjacent normal metal where the dis- 
sipation of spin current provides a non-local mechanism 
to the damping. The damping enhancement is found 
to be proportional to spin mixing conductance, a quan- 
tity playing key roles in the magneto-electronic circuit 
theory.^'* 

The pumped spin current Ip oc M x M is always in the 
plane formed by the free layer magnetization direction M 
and the instantaneous axis about which the magnetiza- 
tion precesses. Therefore, in a single spin valve, when M 
is precessing around the source (drain) magnetization m, 
the pumping current is always in the plane of m and Af.^ 
Let us assume an azimuth angle 9 between m and M. In 
such an in-plane configuration, the pumping current Ip 
has a component Ip sin 9 that is parallel to m. The spin 
transfer torque acting on the source (drain) ferromagnet 
m is the component of spin current that is in the plane 
and perpendicular to m. To simplify the discussion, we 
consider it to be completely absorbed by m. The lon- 
gitudinal (to m) component experiences multiple reflec- 
tion at the source (drain) contact, and cancels the damp- 
ing torque by an amount proportional to Ip sin^ 9 but is 
still aligned along the direction of M x M . Therefore 
the total damping parameter has an angle 9 dependence 
but still picks up a scalar (isotropic) form. This is the 
well-known dynamic stiffness explained by Tserkovnyak 
et al^ In the most general case, when the precessing axis 
of the free layer is mis- aligned with m, there is always an 



out-of-plane pumping torque perpendicular to the plane. 
In the paradigm of Slonczewski, this out-of-plane compo- 
nent is not absorbed at the interface of the source (drain) 
ferromagnetic nodes, while the conservation of angular 
momentum manifests it as a damping enhancement that 
shows the tensor form when installed in the LLG equa- 
tion. 

Studies in lateral spin-flip transistors have suggested 
a tensor form for the enhanced damping parameters 
In spin valves, works based on general scattering the- 
ory have discussed the damping in the framework of 
fiuctuation-dissipation theoremii and shown that the 
Gilbert damping tensor can be expressed using scattering 
matrices,''^ thus enabling first-principle investigation.— 
But explicit analytical expressions of the damping ten- 
sor, its dependence on the magnetic configuration as well 
as the material properties and particularly its impact on 
the magnetization relaxation are largely missing. 

In this paper, we investigate the Gilbert damping pa- 
rameters of the free layer in the so-called dual spin valve 
(DSV).--"-l£ We analyze the origin of the damping tensor 
and derive explicit analytical expressions of its non-local 
dependence on the magnetic configuration and materials 
properties. A generalization of our damping tensor to a 
continuous magnetic texture agrees well with the results 
in earlier works. Particularly, we show, in numeric sim- 
ulations, that by tuning the magnetic configurations of 
the entire DSV, the relaxation time of the free layer can 
be increased or decreased. 
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FIG. 1: A dual spin valve consists of a free layer (with magne- 
tization direction M) sandwiched by two fixed ferromagnetic 
layers (with magnetization directions rriL and mjj) through 
two normal metal spacers. The fixed layer are attached to 



reservoirs. 
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To analyze the spin and charge currents in a DSV, we 
employ the magneto-electronic circuit theory and spin 
pumping,—!^ in combination with diffusion equations.— 
Pillar-shaped metallic spin valves usually consist of 
normal-metal {N) spacers much shorter than its spin-flip 
relaxation length, see for example Ref . [sifisll . To a good 
approximation, in the N nodes, a spatially homogeneous 
spin accumulation is justified and the spin current (J^) 
conservation dictates li = (where subscript i indi- 
cates the source of spin current). 

A charge chemical potential (/x) and a spin accumula- 
tion (s) are assigned to every F or N node. In a transi- 
tion metal ferromagnet, a strong exchange field aligns the 
spin accumulation to the magnetization direction. At ev- 
ery F\N interface, the charge and spin currents on the N 
side are determined by the contact conductance and the 
charge and spin distributions on both sides of the con- 
tact. For example, at the contact between the left lead 
ferromagnet to the left normal metal Ni, called L\Ni 
thereafter, the currents are^ 

Il =^^l [{f-ii - Ml) + - sl) ■ mz,] , 

G 
8n 

+t]l{si - Sl ■ mLrriL)] ■ (1) 

We have used the notation G = + is the sum of 
the spin-cr interface conductance g"^ . The contact polar- 
isation P = (gf" - g^)/{g^ + g^). The ratio 77 = 2g^^/G 
is between the real part of the spin-mixing conductance 
g'^^ and the total conductance G. The imaginary part 
of g"^^ is usually much smaller than its real part, thus 
discardedfi^' The spin-coherence length in a transition 
metal ferromagnet is usually much shorter than the thick- 
ness of the thin film^i^, which renders the mixing trans- 
mission negligible The precession of the free layer mag- 
netization M pumps a spin current Ip — {h/4TT)g'l^M x 
M into the adjacent normal nodes iVi and N2, which 
is given by the mixing conductance g^ at the -F'|A'^i(2) 
interface (normal metals spacers are considered identical 
on both sides of the free layer). 

A back flow spin current at the F|iVi interface reads 

Ii=-^ [2Pf(mi - Mf)M + (si - sf) ■ MM 
Svr 



-?/f(si - Sl • MM)] 



(2) 



on the A^i side. Therefore, a weak spin-flip scattering 
in iVi demands II + Ii + Ip = 0, which is dictated by 
angular momentum conservation. The same conservation 
law rules in N2, where In + I2 + Ip ^ 0. 

For the ferromagnetic (F) nodes made of transition 
metals, the spin diffusion is taken into account properly.^ 
In a strong ferromagnet, any transverse components de- 
cay quickly due to the large exchange field, thus the 
longitudinal spin accumulation s^ = s^rrii, (with v = 
L, R, F) diffuses and decays exponentially at a length 
scale given by spin diffusion length (Xgd) as V^Sj. = 



Sv/^sd- The difference in spin-dependent conductivty 
of majority and minority carriers is taken into account 
by enforcing the continuity of longitudinal spin current 
m^ . ^ -[DlV^sl - DiV^si) at the every F\N in- 
terface. We assume vanishing spin currents at the outer 
interfaces to reservoirs. 

The diffusion equations and current conservation de- 
termine, self-consistently, the spin accumulations and 
spin currents in both TV and F nodes . We are mainly 
concerned with the exchange torque^ T = —M x {Il + 
In) X M acting on M . A general analytical formula is 
attainable but lengthy. In the following, we focus on two 
scenarios that are mostly relevant to the state-of-the-art 
experiments in spin valves and spin pumping: (1) The 
free layer has a strong spin flip (short Xsd) and the thick- 
ness dp > Xsd, for which the permalloy (Py) film is an 
ideal candidate*^ (2) The free layer is a half metal, such 
as Co2MnSi studied in a recent experiment." 

Strong spin flip in free layer. We assume a strong spin 
flip scattering in the free layer i.e., dp > Xsd- We leave 
the diffusivity properties in the lead F nodes arbitrary. 
The total exchange torque is partitioned into two parts: 
An isotropic part that is parallel to the direction of the 
Gilbert damping M x M and an anisotropic part that is 
perpendicular to the plane spanned by mL(R) and M (or 
the projection of M x M to the direction mL[Fi) x M), 
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where the material-dependent parameters "^^f}^'^ and 

San^^ are detailed in the Appendix [Al 

Most interest is in the anisotropic damping described 
by a symmetric tensor with elements 



4«J — 



-m.imj 



(4) 



where i,j — x,y, z (we have omitted the lead index L or 
R). The elements of Aan are given in Cartesian coordi- 
nates of the source-drain magnetization direction. The 
anisotropic damping appears as M x AanM that is al- 
ways perpendicular to the free layer magnetization direc- 
tion, thus keeping the length of M constant. ^1 It is not 
difficult to show that when M is precessing around m, 
the anisotropic part vanishes due to AanM — 0. 

We generalize Eq.Q to a continuous magnetic texture. 
Consider here only one-dimensional spatial dependence 
and the extension to higher dimensions is straightfor- 
ward. The Cartesian component of vector U = M x 
AanM is Ui = —eijkMjTnkmiMi (where Sijk is the Levi- 
Civita tensor and repeated indices are summed). We as- 
sume the fixed layer and the free layer differ in space by 
a lattice constant oq, which allows « Mk{x + oq). A 
Taylor expansion in space leads to U = ~a^M x (SAf), 
where the matrix elements Tiki = {dxM)k(dxM)i and we 
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have assumed that the magnetization direction is always 
perpendicular to d^M. In this case, three vectors d^M, 
M X dxM and M are perpendicular to each other. A 
rotation around M by 7r/2 leaves M and M unchanged 
while interchanging dxM with M x dxM, we have 



^ki = (M X dxM)k{M X dxM)i, 



(5) 



which agrees with the so-called differential damping ten- 
sor Eq.(II) in Ref.;2I]. 

Eq.® suggests that the total exchange torque on the 
free layer is a linear combination of two independent ex- 
change torques arsing from coupling to the left and the 
right F nodes. This form arises due to a strong spin- 
flip scattering in the free layer that suppresses the ex- 
change between two spin accumulations Si and S2 in the 
N nodes. In the pursuit of a concise notation for the 
Gilbert form, the exchange torque can be expressed as 
T = M X with a total damping tensor given by 
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The damping tensor ^ is determined by the entire mag- 
netic configuration of the DSV and particularly by the 
conductance of F|iV contacts and the diffusive proper- 
ties the F nodes. 

Half metallic free layer. This special while experimen- 
tally relevant^S case means Pp = 1- Half-metallicity 
in combination with the charge conservation enforces a 
longitudinal back flow that is determined solely by the 
bias current: The spin accumulations in N nodes do 
not contribute to the spin accumulation inside the free 
layer, thus an independent contribution due to left and 
right leads is foreseen. We summarize the material spe- 
cific parameters in the Appendix [X] When spin fiip is 
weak in the source-drain ferromagnets, « leads to 
"Dis ~ 0. In this configuration, by taking a (parallel or 
anti-parallel) source-drain magnetization direction as the 
precessing axis, the total damping enhancement vanishes, 
which reduces to the scenario of = 1 in Ref . f9| . 

Magnetization relaxation. To appreciate the impact of 
an anisotropic damping tensor on the magnetization re- 
laxation, we perform a simulation, for the free layer mag- 
netization, using Landau-Lifshitz-Gilbert (LLG) equa- 
tion augmented by the tensor damping, i.e.. 
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= -7M X Heff + aoM X 
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M X W . 

dt 
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ao is the (dimensionless) intrinsic Gilbert damping pa- 
rameter. Symbol 7 is the gyromagnetic ratio, Mg is the 
saturation magnetization, and V is the volume of the 
free layer. ^0 stands for the vacuum permeability. The 
dynamics under the bias-driven spin transfer torque is 
not the topic in this paper, but can be included in a 
straightforward wayi^ We give in the Appendix |B] the 
expressions of the bias-driven spin torques. 



We are mostly interested in the relaxation of the mag- 
netization, instead of particular magnetization trajecto- 
ries, in the presence of a tensor damping. The follow- 
ing simulation is performed for the scenario where the 
free layer has a strong spin fiip, i.e.jCase (1). We em- 
ploy the pillar structure from Ref.[l5| while consider- 
ing the free layer (Py) to be 8nm thick (a thicker free 
layer favors a better thermal stability.^^) The source- 
drain ferromagnets are cobalt (Co) and we expect the 
results are valid for a larger range of materials selec- 
tions. The Py film is elliptic with three axes given by 
2a = 90 nm, 2b = 35 nm^ and c = 8 nm. The de- 
magnetizing factors y z in the shape anisotropy en- 
ergy E^,^ = (1/2)^^0^21/ E.=,.y.. ^^Mf are 9x = 0.50, 
Sjy = 0.37 and S>z = 0.13. An external field Ha leads 
to a Zeeman splitting E'zoc = —V^oMgHa ■ M . For Py 
films, we neglect the uniaxial anisotropy. The total free 
energy Ft — Ezcc + -Edcm gives rise to an effective field 
Jfe// = -{l/VM,^iQ)^ET/^M. 

The spin-dependent conductivities in the bulk of Co 
and the spin diffusion length Xco ~ 60 nm are taken from 
the experimental data.^"' For Py, we take Apy « 4 nmj^ 
To have direct connection with experiments, the above 
mentioned bare conductance has to be renormalized by 
the Sharvin conductance,^^ For Py/Cu the mixing con- 
ductance, we take the value g]^S~^ w 15 nm"^,-^ which 
gives D\l{r)f « 1-0. 



(a) B = 50 G; I. P. y-axis. 



(b) Relaxation time 
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FIG. 2: (Color online) Mz as a function of time (in ns) in 
presence of different source-drain magnetic configurations and 
applied fields, (a) The external magnetic field Bz — 50 Gauss 
is applied along z-axis. The blue (dashed), red (solid) and 
black (dotted dash) curves correspond to source-drain magne- 
tization in configurations {y,y), {x,x), and {z,z) respectively, 
(b) Magnetization relaxation times (in the unit of ns)versus 
source-drain magnetic configurations at different applied field 
along 2-axis.: Bz = W G (red □), S;. = 50 G (blue Q), 
Bz = 200 G (green v), Bz = 800 G (black <))• Lines are a 
guide for the eyes. The initial position (LP.) of the free layer 
is taken along y-axis. 

The relaxation time is extracted from the sim- 
ulations by demanding at a specific moment Tr the 
\Mz — 1.0| < 10~^, i.e., reaches the easy axis. In the 
absence of bias, panel (a) of Figl2] shows the late stage 
of magnetization relaxation from an initial position {y- 
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axis) in the presence of an tensor damping, under various 
source-drain (SD) magnetic configurations. The resuhs 
are striking: Under the same field, switching the SD con- 
figurations increases or decreases r^. In panel (b), the 
extracted relaxation times Tr versus SD configurations 
under various fields are shown. At low field = 10 G 
(red □), when switching from {z,z) to {y-,y), Tr is im- 
proved from 8.0 ns to 6.3 ns, about 21%. At a higher 
field Bz — 800 G (black 0), the improvement is larger 
from 5.2 at {z,z) to 3.6 at {y,y), nearly 31%. To a large 
trend, the relaxation time improvement is more signifi- 
cant at higher applied fields. 

In conclusion, combining conservation laws and 
magneto-electronic circuit theory, we have analyzed the 
Gilbert damping tensor of the free layer in a dual spin 
valve. Analytical results of the damping tensor as func- 
tions of the entire magnetic configuration and material 
properties are obtained. Numerical simulations on LLG 
equation augmented by the tensor damping reveal a tun- 
able magnetization relaxation time by a strategic selec- 
tion of source-drain magnetization configurations. Re- 
sults presented in this paper open a new venue to the 
design and control of magnetization dynamics in spin- 
tronic applications. 

X.Wang is indebted to G. E. W. Bauer, who has 
brought the problem to his attention and offered invalu- 
able comments. 



spin flip in free layer. 
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Appendix B: Bias dependent spin torques 



Appendix A: Material dependent parameters 

In this paper, ^l(r)f = .9l{_r)/5f^ mixing con- 

ductance ratio and Xl(r) = ''^l(r) ' The diffusivity 
parameter ^^.(fl) = (I>l(r){'^~ Phm) hL{R), where for the 



left F node 



1 



1 



ihSal cr'l tanh(di / A ^ ) 



Gl{1-PI) 



(Al) 



where h the Planck constant, S the area of the thin film, 
e the elementary charge, Al the spin diffusion length, c?l 
the thickness of the film, and cr^'-''' the spin-dependent 
conductivity. 0^ is obtained by substituting all L by i? in 
Ea. (|Aip . Parameter is given by S^p = [1 — Pp)(j)p/T]p 
with 
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(o^+a^O^Ffl^^^i p2^ 
-L + ,.o„T„4 Lxf(1 - ^p) 
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The material dependent parameters as appearing in the 
damping tensor Eq.® are: (1) In the case of a strong 



The full analytical expression of bias dependent spin 
torques are rather lengthy. We give here the expres- 
sions, under a bias current /, for symmetric SD fer- 
romagnets (i.e., 4>l = 4>R ^ 4> thus = ^fl. = ^) 
with parallel or anti-parallel magnetization direction. (1) 
With a strong spin fiip in the free layer, the parallel 
SD magnetization leads to vanishing bias-driven torque 
T"!^' = 0; When the SD magnetizations are anti-parallelly 
(i.e., niL = -ruR = m), 

-f-<n(fn + (,f{1- x^) + X^)] mpx{mx mp). 

(Bl) 

(2) When the free layer is half metallic, for symmetric 
SD ferromagnets , Tj^'' — and 

rrib) Ih (j)P 

J-AH = T TTT ^TT T7-^ XTT'TIp X [171 X mp). 

(B2) 
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